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Abstract: We present a theory of vortex liquid-to-solid transition in homogeneous quasi 
2D superconductors. The free energy is written as a functional of density of zeroes of the 
fluctuating order parameter. The transition is weakly first-order and well below the H C 2(T) 
line. Transition temperature, discontinuities of the average Abrikosov ratio and of the average 
superfluid density, the Debye- Waller factor and the latent heat are in good agreement with 
Monte Carlo simulations. The density is only weakly modulated in the "vortex-solid" phase, 
consistent with the density-wave behavior. 
PACS: 72.20.My, 71.45.-d, 71.30.+h 



The experimental observation of the wide region around the H&iT) line where the 
Abrikosov lattice appears melted demonstrates the importance of thermal fluctuations in 
high temperature superconductors (HTS). Such fluctuations are generally important in quasi 
two-dimensional (2D) materials, such as thin films or layered systems with weak interlayer 
coupling. The existence and the nature of ordered low-temperature phase is still matter of 
some controversy. The early theories [Q] assumed the ordered phase at low temperatures and 
proposed dislocation-unbinding as a mechanism of second-order melting transition or used 
the Lindemman criterion to locate the melting point ||. However, even the very existence 
of vortex-lattice phase, which is the starting point of these theories, cannot be taken for 
granted @, ||. External magnetic field causes the system to appear D-2 dimensional on 
perturbative level, and the perturbation theory carried out to eleventh order |J shows no 
sign of low-temperature Abrikosov phase in 2D. Renormalization group analysis near the 
upper critical dimension D=6 found the continuum of relevant charges to be generated by 
the renormalization procedure and indicated that the transition into Abrikosov vortex-lattice 
state might be first order, as expected from the simple Landau argument 0. The relevance 
of this analysis for real superconductors, however, is not obvious since the lower critical 
dimension for the Abrikosov transition is D=4. Recently, Monte Carlo (MC) simulations 
|8|, |9], [Hj revealed very weak first order vortex liquid to solid transition in 2D. At present, 
our information about this transition comes almost exclusively from such numerical work. 

The main purpose of this paper is to develop a working analytic theory which has 
predictive power and which can be applied to a wide variety of experimentally relevant 
issues. We derive a density-functional representation of the free energy for zeroes of the 
fluctuating superconducting order parameter, with the exact two-body correlations built in. 
Our theory is a "first principle" in the sense that it can be obtained from the microscopic 
description by a series of well understood approximations. It is then shown that above 
certain values of the first two peaks in the vortex structure factor triangular lattice has lower 
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free energy than the uniform configuration. Using the information on the dependence of the 
structure factor on temperature from the existing MC simulations we locate the transition 
temperature and calculate several typical quantities like the latent heat and the change in 
thermal average of the Abrikosov's ratio at the transition. Obtained results are in good 
agreement with those found in numerical simulations. 

The starting point is the Ginzburg-Landau (GL) partition function for a homogeneous 
2D superconductor in perpendicular magnetic field, with fluctuations of the magnetic field 
neglected (k ^> 1). We are interested in the regime where the Landau level (LL) structure 
of Cooper pairs dominates the fluctuation spectrum: This is the case for fields above « 
{6 1 16) H c2(0) (T/Tcq), where 9 is the Ginzburg fluctuation parameter fll] (For example, in 



BSCCO 2:2:1:2, 9 ~ 0.045 and ~ 1 Tesla). In this regime, the approximation in which 
one retains only the lowest Landau level (LLL) modes describes essential features of the 
problem. The partition function is then Z = J D\tp*ip] exp (-S), and 

8 = i{a\T) [ d 2 r\i,{f)\ 2 + | f d 2 f\m\ 4 } , (1) 
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where ot'(T) = a(T)(l — H/H^T)), d is the thickness of the film or the effective interlayer 
separation and at(T) and (3 are phenomenological parameters. We choose to work in the 
symmetric gauge so that the order parameter is a holomorphic function and can be written 
as ^Jj(z) = 4>Y[f = i(z — Zi) exp — (|z| 2 /4) where N is the area of the system in units of 2tt1 2 , 
z = (x + iy)/l and / is the magnetic length for the charge 2e. The partition function can now 
be expressed in terms of variables <fi and {z^} ||. They represent two distinct tendencies 
in a superconductor; describes the overall growth of the local superconducting order, 
while {z^ represent the remaining weak lateral correlations between vortices. If the latter 
are treated in a certain average way one can account for 98% of the thermodynamics [[IT 



However, the transition that we want to discuss is entirely within the remaining 2%, and the 
correlations among {z{\ are crucial. The integration over <ft can be performed exactly in the 
thermodynamic limit N — » oo, yielding Z = j JJ^dzidz* /2n)f i N2 Jl\zi — Zj\ 2 exp{— S'}, 



where 

S ' = ~T + \ {V2 + 2)1/2 + sinh_1 ( F /^))> ( 2 ) 
V({zi}) = gP/yfP, P({zi}) = J(dzdz*/2rcN) U* \z- Zl \ n exp ~(n\z\ 2 /4) and g = a'^l 2 d/(3T. 
The variable <fi is determined by 

< H 2 > P 2 ^^ = v\i + VT+W^). (3) 

The original problem is now equivalent to the thermodynamics of classical 2D system of 
particles which we call dense- vortex-plasma (DVP) ||. We note several features of DVP 
system which should be reflected in our density-functional: 1) particles interact via long- 
range multiply-body forces, 2) the system is incompressible, 3) the system is scale invariant, 
4) the thermodynamics depends on a single dimensionless coupling constant g. Previous MC 
simulations show transition from a liquid to a solid state taking place at gF ~ —7 || [|, |10j. 
The mechanism of the transition, however, remains obscure in the partition function (2) 
since the exponentiated energy is still a perfectly smooth function of V({zi}). From Eq. (2) 
we can write the exact free energy per vortex as f(g)/T = S' — s(< (3a >) where < (3a > is 
the ^-dependent thermal average of the Abrikosov ratio (3a{{zi}) = / 4 /(/ 2 ) 2 determined by 
minimizing f(g). The entropy s(< (3a >) is given by 

S (< (3 A » = at- in/ n d -^fr N/2 n i* - %i 2 ^- < & »• (4) 

i i<j / 

The transition is (as usual) hidden in the entropy part of the free energy, but the evaluation 
of the integral (4) is a hopeless task. Obviously, a different route needs to be taken. 

To study solidification transition it is often beneficial to change variables from parti- 
cle coordinates to density of particles. This can be achieved by inserting the unity 1 = 
/ D p(r)5[p(r) — Z)5(r — r*j)] in Z (2). After the 5-function is expressed as an integral over 
auxiliary field $(r) coordinates of vortices can be integrated out. The partition function 
becomes Z = J Dp{r)D<&{r) exp{— 5"'} , 

S" = U[p(f)} - J d 2 r$(r)p(r) -N\nJ d 2 fexp (-$(f)) (5) 
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where U[p] is the energy density-functional determined by the details of the DVP interaction 
and we changed variables z$ — > $. We expand U[p] in terms of multy-body interactions and 
keep only the first, two-body term so that U = ~ / p(f[)v{f\ — r^)p(r^) + 0(p 3 ). After 
this truncation of the energy functional, density can be integrated out leaving the partition 
function as an integral over field $(r) with the action 

F = — - J dr\dr 2 §{fi)v~ l (f{ - r^)$(r^) - N\n J d 2 rexp (-$(f)). (6) 

We take the Eq. (6) to be the desired expression for the mean-field free energy of the vortex 
system. In order to account for the strong two-particle correlations that exist in the liquid 
state of DVP we make the substitution —v(r) = c(r), where c(r) is the direct (Ornstein- 
Zernike) correlation function of the vortex-liquid whose Fourier transform is related to the 
structure factor via the relation S(q) = (1 — c(g)) -1 [12|. The structure factor S(q) is the 
Fourier transform of the connected density- density correlation function. 

Before proceeding with the analysis of the phases described by the free energy (6) let 
us comment on validity of the truncation of the DVP energy density functional. Even if the 
particles interacted only via two-body forces all multi-body correlations will be generated 
and the true free energy would contain additional terms in (6) of higher order in field 
Therefore, the fact that particles in DVP interact through multi-body forces makes the 
approximate free energy (6) no less appropriate here than for the ordinary correlated liquid. 
In fact one might expect that the long-range interaction operating in DVP will make the 
system only more mean-field like. The important feature of the approximate free energy (6) 
is the last term which is purely entropic and hence crucial for the freezing transition. The 
free energy (6) may be thought of as a harmonic approximation around the liquid state and 
it will be shown that this particular form of it is very good for an incompressible system. 

We assume that the freezing of the DVP happens into triangular lattice with period 
a, a 2 = An / \/3 and we work in units I = 1. Note that this assumption is not necessary 
and that we could consider other lattices. Actually, at this point it is not obvious that 



the ordering wave-vector will be in any simple relationship to the magnetic length. Within 
the resolution of the MC data for the structure factor, however, the strongest tendency 
for ordering corresponds to the triangular modulation. We determine the solid phase by 
two real order parameters $1 and $ 2 in the following way: $(f) = $ — 2 Y%=i ®ifi(r)i 
where fi(r) = 1/2J2gi expiGif = cos(arc) + 2 cos(a-\/3y/2) cos(aa;/2) and the sum runs over 
six shortest reciprocal lattice vectors \G\ = G\ = a. The same sum, but over next six 
shortest reciprocal vectors \G\ = G 2 = V3a determines the function f?,(r) = cos(ay/3y) + 
2 cos(3ax/2) cos(av / 3y/2). Inserting this into expression (6) we get two-order-parameter 
form of the free energy per vortex 

F = |i + $ + 3^^-ln/ d 2 f exp (2 ]T $^(f)) (7) 
2c i=1 a J i=1 

where ci/ 2 = c(Gi/ 2 ). The free energy is minimized by the solutions of the equations: 

3$, = /d 2 f > / J -(f)exp(2S? = i^/»(^) (g) 
Cj /«Prexp(2^ =1 $ i / i (f)) {) 

for j = 1, 2. The equation for the uniform component is $o = ~ Co and system is manifestly 
incompressible since the average density p(G = 0) = 1 does not change with temperature. 
It is easily seen that / fj(r) = and one solution of the equations is always $i = $ 2 = 
which describes the liquid state of DVP. To find a non-zero solution, we choose a pair 
(ci,c 2 ), numerically solve Eqs. (8) and compare free energies of the solid and the liquid 
configurations. This way one obtains the crosses joined joined by the full line with a negative 
slope on the structural phase diagram presented on Fig. 1. For values of Si and S 2 above the 
line triangular lattice of vortices is the stable phase, while below it DVP is in liquid phase. 

To find the value of the coupling constant gp where the DVP freezes the information 
on Si and S 2 as functions of g is needed. The vortex structure factor can be found by using 
the perturbation theory to determine the structure factor for the superfluid density, |?/>(r)| 2 , 
and then connecting the two through the expression for tp(r) given below Eq. (1) and some 
of the standard correlation function hierarchies (BBGKY, hypernetted chain, etc.). While 



this is possible in principle it is more expedient here to utilize the vortex structure factors 
from MC simulations || 13|]. Particularly useful is a detailed analysis of the peaks in the 
vortex structure factor by O'Neill and Moore Ul3| . In the MC simulations in the spherical 
geometry they find that, for 5 < g 2 < 50, both peaks depend approximately linearly on g 2 . 
The straight line with positive slope in Fig. 1 represents Si — 5*2 relation derived from there. 
Two lines intersect at S[ = 4.45 and S£ = 1.47, which gives gp = —6.5 and the values 
of order parameters at the transition are $i = 0.50 and $2 — 0.10. Interestingly, around 
this value of the coupling constant O'Neill and Moore first start to see critical slowing down 
of the dynamics in their MC simulations, but they argued against finite temperature phase 
transition. Our result is in excellent agreement with MC simulations in Refs. f| |9|, |10| . 

The main features of DVP (long-range interactions, incompressibility) resemble those 



of another well studied classical system, namely 2D one component Coulomb plasma [14 



It is therefore not surprising that the values of first two peaks of the structure factor that 
we found at the transition closely match the numbers found in MC simulations there. In 
particular, the value of S[ is very close to 4.4 which is known as 2D version of the Verlet 
criterion. For comparison, our number for is much worse for the hard-disk system where 
it equals 1.9 ||14[| . This is because the system of hard-disks is compressible and the proposed 
form of the mean-field free energy (6) ignores the change in average density at the freezing 
transition. 

Since the transition takes place at a finite value of Si the transition is first order in 
agreement with what has been observed in recent experiments in clean YBaCuO single 
crystals [fl5 1 . The latent heat comes solely from the structural change at the transition. 
If we take the expression for entropy of the system to be s = J d 2 fp(f)ln[p(f)} the latent 
heat equals Q = NTpS^l/cf + $2/°^)- For the above numbers this gives Q/NTp = 1.0, 



somewhat higher than the number 0.4 found in MC simulations P, 10, [L4|. Going back to 



the exact expression for the free energy it follows that the thermal average of Abrikosov ratio 
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has a discontinuity at the transition A < (3a > / < (3a >~< (3a > AQ/g F NT F = 0.025 if 
we take < (3a >~ 1.2 at g — g F . Invoking Eq. (3) we see that the average superfTuid density 
has a discontinuous jump of about 2%. This in reasonable agreement with the results of 
Refs. |§ and [10|. The Fourier components of the density of vortices can be determined 
from : 

-> Jd 2 fexp{iGf - $(f)) 
P[ ' /d 2 fexp(-$(f)) { ' 

and the Debye- Waller factor is given by u(G) = |p(G)| 2 . We find that for \G\ > G\ p(G) « 

0.72 exp(— X 2 G 2 ) with A = 0.47. Similar Gaussian fall-off was found for freezing into 3d bcc 

or fee lattices, but with smaller coefficient A (\b cc = 0.34, A/ cc = 0.19, jT2"| ). The fast decay 

of Debye- Waller factor indicates that the transition is very weakly first order and the density 

modulation is rather weak right below the transition. 

The present scheme allows us to examine the possibility of high temperature solid phase. 
We determined the values of (c\,C2) where the local $i ^ 0, $2 7^ minimum of the free 
energy first appears; at that point free energy of the liquid state is still lower than the one of 
the solid. This way triangles joined by the dashed line on Fig. 1 are obtained. Intersection 
with the S\ — S2 line then determines the highest value of g where the solid is still locally 
stable: g s h = —6.25. The interval Ag sh = g sh — g F = 0.25 agrees well with the value 0.2 
obtained in Ref. |§ for the finite size system. The supercooling of the liquid phase is also 
possible, and the lowest g = g sc where the liquid phase is locally stable is given by the value 
where the first peak of S(q) diverges. From the data in Ref. || we expect that the interval 
Ag sc = g F — g sc is of the same size as Ag sh . The supercooling and the superheating data 
from the same reference indicate that there is a well defined critical value of Abrikosov ratio 
< (3a >~ 1.18 at which the DVP system can undergo a continuous transition. 

The vortex-liquid freezing transition discussed here is not a superconducting transition 
in the sense of breaking of U(l) symmetry ||. It is best thought of as a transition to a 
charge- density- wave of Cooper pairs. This is consistent with the weak density modulation of 
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the ordered phase. The superconducting pairing susceptibility remains short-ranged in this 
density-wave, although this range becomes 3> I as density modulation increases. This phase 
should be distinguished from other phases described in the literature [|16|]. True supercon- 
ducting transition in high magnetic field is almost always due to disorder, and, unless the 
disorder is very strong, it will be much below the freezing transition discussed here. Finally, 
the weak interlayer coupling in HTS and similar systems will act to "lock-in" density-waves 
from different layers resulting in a 3D-ordered weakly modulated structure. The first order 
character of the transition will be enhanced. The change in the transition temperature will 
be rather small and can be easily estimated ||. 

To conclude, we have presented a density-functional theory of the vortex liquid-to-solid 
phase transition in homogeneous quasi 2D superconductors. As the input, we have used 
the temperature-dependent first two peaks of the numerically determined structure factor 
for the zeroes of the superconducting order parameter. As the output, we obtain the value 
of transition temperature, the discontinuities of thermal average of the Abrikosov ratio and 
of the average superfluid density, the Debye- Waller factor at the point of transition and 
the latent heat. We determine the temperature interval for superheating of the vortex-solid 
phase and conjectured that there is a critical point that can be reached by supercooling. 
The results are found to agree well with Monte Carlo simulations. While we have used 
the numerical results as an input to our theory, it is not necessary to do so. Large order 
perturbation expansion results for the structure factor of |\l/(r)| 2 [[L7|] could have been used 
instead so that the theory is kept entirely at the analytic level. 

This work has been supported in part by the David and Lucile Packard Foundation. 
One of us (IFH) thanks Professor M. O. Robbins for useful discussions in the early stage of 
this work. 
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Captions: 

Figure 1. Structural phase diagram for vortex system in terms of first two peaks of 
the structure factor of vortices. Crosses denote the values of Si and S 2 where the vortex- 
solid becomes energetically favorable, and triangles mark the points where the solid first 
becomes locally stable. The line with positive slope depicts Si — S2 relation inferred from 
MC simulations [13j . The transition occurs at the intersection of this line with the line which 
joins the crosses. 
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